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Abstract 

4C +21.35 is a flat-spectrum-radio-quasar-type blazar, in which the rapid vari¬ 
ability of very high energy (VHE, ^ 100 GeV) emission as short as ~ 10 minutes 
was observed by MAGIC Cherenkov telescopes, and the VHE spectrum extends up 
to at least 400 GeV. In this paper, by using a flat broad-line region (BLR) structure, 
we study the location and properties of y-ray emitting region of 4C +21.35 under the 
constraints of multiwavelength data. We £t three quasi-simultaneous spectral energy 
distributions (SEDs) using homogeneous one-zone leptonic model, in which the flat 
BLR with the aperture angle of a = 25° and a spherically symmetric hot dusty torus 
with the temperature of Tsub = 1200 K, are assumed. The results show that the jet 
structure of 4C +21.35 is nearly conical with a half-opening angle of 6j ~ 0.29° — 0.6°. 
Furthermore, the emitting region is located within the BLR clouds and approaches 
to outer radius of the BLR during the flaring states, while it is well beyond the 
dusty torus in quiescent state. The quiescent high-energy emission is dominated by 
synchrotron self-Compton (SSC) process, the high-energy emission during the flaring 
periods is dominated by Compton scattering of BLR and dusty torus photons with 
the value of Compton-dominance parameter is about 30. Moreover, the fit to opti¬ 
cal/ultraviolet data provides a further support that the central black hole (BH) mass 
of 4C +21.35 is 6 X 10® M©. 

Key words: galaxies: active - gamma rays: general - quasars: individual (4C 
+21.35) - radiation mechanism: non-thermal 
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1. Introduction 


Only ~ 10 — 20% radio-loud active galactic nuclei (AGNs) have jet close to our line of 
sight, i.e. blazars (Ghisellini 2011). Their emissions are dominated by non-thermal emission 
throughout entire electromagnetic spectrum. Blazars are an outstanding subclass of AGNs, 
which consists of BL Lacertae objects (BL Lacs) and flat spectrum radio quasars (FSRQs) 
(Urry & Padovani 1995). 

BL Lacs and FSRQs have signihcant differences in their structures. BL Lacs could 
have a radiatively inefficient accretion disk and BLR (Stocke et ah 2011), but lack of dusty 
torus. Instead, FSRQs have a radiatively efficient accretion disk, together with BLR and dusty 
torus. The emission from BL Lacs is generally dominated by synchrotron self-Gompton (SSG) 
process (Finke et ah 2008; Dermer et al. 2009), while the emission from FSRQs is dominated by 
the external radiation Gompton (ERG) mechanism, in which the seed photons could originate 
from the accretion disk (Dermer et al. 1992; Dermer & Schlickeiser 1993), BLR (Sikora et al. 
1994; Donea & Protheroe 2003; Liu & Bai 2006; Gelotti et al. 2007; Ghisellini & Tavecchio 
2008; Ghisellini & Tavecchio 2009; Tavecchio & Ghisellini 2012; Lei & Wang 2014) and dusty 
torus (Donea & Protheroe 2003; Blazejowski et al. 2000) radiation helds. 

4G -1-21.35 (PKS 1222-1-216, z=0.432; Osterbrock & Pogge 1987) is a y-ray-emitting 
FSRQ (Hartman et al. 1999), conhrmed by the Energetic Gamma Ray Experiment Telescopes 
(EGRET). During the hrst three months of scientihc operation of the Fermi Large Area 
Telescope (Fermi-LAT), the flux of this source was too low to be included in the Fermi-LAT 
Bright Source List (Abdo et al. 2009). Fortunately, it was included in the First LAT Gatalog 
(Abdo et al. 2010), which was associated with IFGL J1224.7-1-2121. Fermi-LAT observations 
show that 4G -1-21.35 was in a quiescent state during 2008 August to 2009 September and 
subsequently the source became active. On 2010 April 24, the Fermi-LAT detected a strong 
outburst from this source (Donato 2010). During this outburst, the observed y-ray emission 
extends up to the photon energies larger than 100 GeV (Neronov et al. 2010; Neronov et al. 
2011). In 2010 April and June, the source arose a strong GeV outburst, several major flares 
are described by rise and decay timescales of order of day (Tanaka et al. 2011). On 2010 June 
17, the Major Atmospheric Gamma Imaging Gherenkov Telescope (MAGIG) detected y-ray 
emission from 4G -1-21.35, where the flux-doubling timescale of the y-ray was about 10 minutes, 
and the spectrum extends from about 70 GeV up to at least 400 GeV (Aleksic et al. 2011). 

Gurrently, several explanations on the production of the short y-ray variability of 4G 
-1-21.35 have been put forward. Kushwaha et al. (2014a) considered the emitting region to 
be at the recollimation zone of the jet. Tavecchio et al. (2011) proposed three scenarios to 
explain the SEDs and favored the emitting region to be re-collimated and focussed zone of 
the jet at larger distance from the black hole (BH). In addition, photons/axion transition was 
proposed to explain fast y-ray variability (Tavecchio et al. 2012; Galanti et al. 2014). The 
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comprehensive SED investigations to 4C +21.35 have been performed by Ackermann et ah 
(2014) (hereafter, Paper-I), in which three quasi-simultaneous SEDs from 2009 to 2011 were 
modeled with a combination of SSC and Compton-scattered soft photons from dusty torus. In 
addition, Dermer et al. (2012) argued that the rapid VHE radiation originates from a beam of 
leptonic secondaries initiated by the ultra-high energy cosmic rays (UHECRs) at the pc-scale 
region. 

For 4C +21.35, the very rapid y-ray variability indicates an extremely compact emitting 
region, with size < cfvar^D/(l +- 2 ) — 1-26 x 10^®(5D/100)(tvar/10minutes) cm, where Aar is 
the observed flux-doubling timescale and (5 d is the Doppler factor. Here, we have assumed 
that the emitting region has a higher Doppler factor of 100. Under conical conhguration, the 
half-opening angle of relativistic jet is generally assumed to be ~ 0 . 1 , this corresponds to the 
emitting region at radial distance r ~ 10^® cm. Generally, the BLR has the scale of ~ 0.1 — 1 pc, 
this implies that the emitting region is well within the BLR cavity. Under this situation, the 
low-energy photons originating from BLR are the target source for ERC process. Meanwhile, 
they will also severely absorb the 7 -rays above a few tens of GeV via pair production due 
to 7 — 7 interactions (Liu & Bai 2006). The use of a flat BLR will effectively alleviate this 
absorption effect, as presented by Lei & Wang (2014). 

In Section 2, we give a brief description of the model. In Section 3, we perform spectral 
modeling to three SED data points. Our discussion is given in Section 4. The conclusion is 
given in Section 5. Symbols with a numerical subscript should be read as a dimensionless 
number = X/(10’^cgs units). Throughout the paper, unprimed quantities are measured in 
the observer’s frame, unless stated otherwise, and primed quantities refer to the blob’s frame, 
while the starred quantities refer to the stationary frame with respect to the BH. We adopt the 
cosmological parameters: Hq = 71 km s“^ Mpc“^, Dm = 0.3, Da = 0.7. 

2. Model description 

The model is set up within the framework of homogeneous one-zone leptonic model, 
in which for simplicity the emitting region is assumed to be a spherical blob with the bulk 
Lorentz factor Tj. The jet makes an angle 6^ with respect to the line of sight. The blob is 
hlled by tangled, homogeneous magnetic held of strength B, together with relativistic electrons 
(including electrons and positrons). The relativistic effects are characterized by Doppler factor 
(fn = [rj(l —/3^cos0v)]~^ where (3r is the bulk velocity in units of the speed of light c. The 
comoving size of the blob is related to the observed variability timescale tvar through ~ 
AarC5D/(l + ^)- 

A broken power-law function with a cutoff is adopted to describe the relativistic electron 
distribution within the emitting region: 



where no is the normalization factor; 7 ' = E'^/rrieC^ is the electron Lorentz factor assnmed to 
vary within the range 7 ^,^ < 7 ' < 7 max- -Si S 2 are the spectral indices below and above the 
break energy 7 (,, respectively. 

An important parameter to characterize 4C +21.35 is the central BH mass, which is still 
a matter of debate. Wang et ah (2004) snggested a valne of 1.48 x 10® Mq obtained by the 
relation between H/3 broad line and the continnnm Inminosity-BLR radius. Recently, based on 
the virial method from the full width at half-maximum of Mgn, H/3 and Ha broad lines and 
thermal continuum luminosity. Farina et ah (2012) presented a larger value of ~ 6 x 10® Mq, 
which agrees well with the values found by Shen et ah (2011) and Shaw et ah (2012). Their 
results originate from Sloan Digital Sky Survey (SDSS) spectra. Given these results, the BH 
mass of 4C +21.35 Mbh = 6 x 10® Mq is adopted in our model. Beyond the BH, there is 
an optically thick accretion disk; its spectral luminosity is assumed to be characterized by a 
Shakura-Sunyaev disk (Shakura & Sunyaev 1973) spectrum in the form (Dermer et ah 2014) 


+^d(e*) 


1.12Ld 
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where e* = E^./rrieC^ is the dimensionless photon energy in stationary frame, which is related to 
observed dimensionless photon energy by e* = (1 + z)e. The value of Cmax is the function of the 
spin of the BH and relative Eddington luminosity. 

Generally, the BLR is assumed to be spherical, and with a radius that can be written in 
terms of the accretion disk luminosity as (Ghisellini & Tavecchio 2008) 


Rblr — 0.1^ 


■pc. 


( 3 ) 


y 10 ^®erg s“^ 

In this paper, we consider a flat BLR, in which the cloud matter is hlled within a region with 
width extending from rBLR,in = (Rblr - 0.05) pc to rBLR,out = (Rblr + 0.05) pc. A detailed 
description of the flat BLR structure is given by Lei & Wang (2014). 

Moreover, the spectral luminosity of the infrared (IR) component can be approximated 
by a graybody radiation held, namely, 

15Lir (e*/0)" 


G-hiR(e*) — 


M exp(e*/ 0 )-l’ 


( 4 ) 


where 0 = k-QT^nh/'me(? ■, the sublimation temperature of the dusty torus of 4G +21.35 is Tgub = 
1200 K (Malmrose et al. 2011). For simplicity, the dusty torus is assumed to have a geometrically 
thin and spherically symmetric structure with radius that corresponds to dust sublimation 
(Malmrose et al. 2011), i.e.. 




torus 


~ 0.4 


/ 

f \ 

1 10 ^® erg s“®' 

U500Ky 


- 2.6 


pc. 


( 5 ) 


Gombined the H/3 luminosity, Lh^ ~ 2 x 10^® erg s and the scaling relation Lblr — 
25.3 Lh^ (Wang et al. 2004; Fan et al. 2006), Tanaka et al. (2011) derived the accretion-disk 
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luminosity as — 10 Lblr — 5 x erg s“^, where the fraction 0.1 of the disk radiation 
reprocessed by BLR clouds has been assumed. Tavecchio et al. (2011) assumed that the Swift 
ultraviolet spectrum of 4C +21.35 is contributed by the accretion-disk radiation, and estimated 
that Ld — 5 X 10^® erg s“^. Moreover, Ackermann et al. (2014) performed the SED htting to 
4C +21.35, and restricted the accretion-disk luminosity to Ld = 0.2 LEdd = 1-56 X 10^® erg s 
We use the larger disk luminosity here, noting that a larger Ld corresponds to a larger BLR 
radius, subsequently reducing the absorption efficiency of the y-rays. 

3. VHE emission production and the SED modeling 

3.1. Constraint on the site of the VHE -ray production 

According to the aforementioned discussions, the y-ray emitting region associated with 
~ 10 minutes VHE variability is roughly located at ~ 1.3 x 10^® cm ~ 5 x 10“^ pc, assuming that 
this region takes up entire cross-section of the jet. In a spherical BLR scenario, the y-rays with 
energies above tens of GeV are expected to be strongly attenuated by the interactions with 
low-energy photons from the BLR. Therefore, it is useful to probe the location of the y-ray 
emitting region based on a flat BLR structure. 

Figure 1 illustrates a calculation of the pair-production opacity for different y-ray ener¬ 
gies from = 50 GeV up to E-^ = 400 GeV, i.e., throughout the whole electromagnetic spectrum 
available by MAGIG from 4G +21.35. Aleksic et ah (2011) have shown that the de-absorbed 
VHE spectrum smoothly connects to the Fermi-LAT spectrum, suggesting that the same pop¬ 
ulation is powering both emissions. Generally, the short variability is thought to be related 
to the compact emission region, and correspondingly with shorter distance from the BH. In 
order to explore the location of fast VHE variability observed on 2010 June 17, the location of 
the emitting region is hxed at rt = 5 x 10“^ pc, corresponding roughly to the short timescale 
variability of 10 minutes. However, it should be emphasized here that in comparison to the 
observation on 2010 June 17, Fermi-LAT observations during its active period showed that the 
source generated a few y-rays flares in the MeV-GeV band with longer variability timescales 
(Foschini et al. 2011; Kushwaha et al. 2014b). 

We can see from Figure 1 that the observed VHE photons can not be produced within 
the BLR cavity. This constraint is also effective for lower y-rays with energy E-^ = 50 GeV, 
where the opacity r remains larger than unity, even for a ~ 10°. Therefore, if the jet of 4G 
+21.35 has a half-opening angle ~ 0.1, it is not appropriate to make the assumption that the 
location of emitting region can be determined by combining the measured variability timescale, 
~ cJoGar, and the assumption that the emitting region covers the entire cross-section of 
the jet. Under this case, the location of the emitting region related to ~ 10 minutes variability 
must be located within the BLR cavity, in direct conflict with the theoretical expectation. This 
seemly indicates that the jet structure of 4G +21.35 is exceptionally complicated, the ultra- 
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short variability could originate from the spine of the structured jet, or from the re-collimated 
and focussed zone of the jet. 

In our scheme, the pair-production opacity r depends on the 7 -ray energy and on two 
geometrical parameters: the BLR aperture angle a and the location of the emitting region r. 
Letting r = 1, we can hnd out the critical position Vc of the emitting region, as shown in Figure 
2. The shaded region represents the part surrounded by inner radius rBLR,in and outer radius 
^"BLR,out of the BLR, whereas the horizonal thick line represents the scale of the dusty torus. In 
calculation, the opacity r is a path integration along the line of sight from Vc to rmax, where 
^max = 5 pc is hxed. It is clear that when the emitting region is within the BLR cavity, just for 
a less than 35°, the 7 -rays with energy = 100 GeV can escape the absorption of the BLR 
photons. For the observed higher-energy 7 -ray photons, the location of the emitting region 
must be within or beyond the BLR clouds, this provides strong constraint on the location of 
7 -ray-emitting zones for FSRQs. 

3.2. Modeling the SEDs of 4C +21.35 

Three quasi-simultaneous SEDs shown in Figure 3, were presented originally by Paper- 
1. The MAGIC data have been corrected for EBL absorption using the model of Finke et al. 
(2010). These SEDs include a quiescent state, in which the temporal flux is integrated over 
time spanning from 2008 August 4 to September 12 (dark-gray diamonds, SED-1), together 
with the flaring states of 2010 April 29 (blue squares: SED-2) and 2010 June 17 (red circles, 
SED-3). 

We model the SEDs of three epochs within the framework of one-zone leptonic scenario. 
At hrst, we obtain information on the peaks of the spectrum and locations at large. In syn¬ 
chrotron hump, three peaks could occur, including the hrst peak at radio band, the second one 
at infrared band from the reprocessed emission of the dusty torus, and the third one at ultravi¬ 
olet band from the emission of the accretion disk. In the second hump, we can not identify the 
number of peak and the corresponding location explicitly. However, two peaks could appear at 
MeV and GeV bands, in which the hrst peak could from the SSG emission and/or the Gompton 
scattering of low-energy photons from accretion disk depending on the distance from the BH, 
and the second peak could originate from the Gompton scattering of the dihuse photons from 
the BLR and/or the dusty torus. In the SED modeling, we due with the radio data as an 
upper limits, because they could be the superposition of multiple self-absorbed components or 
originate from the extended regions of relativistic jet. 

In order to reduce the number of free parameters used in the model, several parameters 
have been hxed based on the following arguments: 

(1) The sizes of the emitting regions for three multiwavelength campaigns are equally 
hxed at = 1 X 10^^ cm, which corresponds to short variability with timescale ~ 10 minutes. 

( 2 ) The measurement of the apparent superluminal motion to relativistic jet of 4G -1-21.35 
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gives a bulk Lorentz factor Fj 20 (Jorstad et al. 2001) on parsec scales, and constrains the 
opening angle within ~ 3°. Moreover, blazar jets have the opening angles that satisfy the 
relation 6^ ~ (0.1 — 0.2)/rj (Pushkarev et al. 2009; Clausen-Brown et al. 2013) on average. Here 
we use a larger jet half-opening angle of 6j ~ 0.2/rj, limiting the jet opening angle of 4C -1-21.35 
to be 6j ~ 0.57°. In SED modeling, we take the value 6j = 0.6° approximately, and the viewing 
angle is eqnal to the jet half-opening angle, 9^ = Oy 

(3) The typical BLR thickness ratio was fonnd to be H/R ~ 0.4 (Decarli et al. 2008; 
Decarli et al. 2011; Shaw et al. 2012), corresponding to a ~ 25° approximately. Thns we hx 
the aperture angle a of the BLR to be 25°. We also assume that the BLR clouds intercept 
a fraction /blr = 0.1 of the central continuum into the diffuse radiation and the dusty torus 
intercepts a fraction /ir = 0.3 of the one into the infrared. Again, we further assume that the 
BLR reprocesses a fraction of /cont = 0.05 and /une = 0.05 into continnnm and emission lines, 
respectively. 

We ht the nltraviolet data with a blackbody spectrum from the standard thin disk, in 
which the accretion-disk Inminosity is Ld = 4.6 x 10^® erg s“^. Indeed, the infrared radiation of 
dusty torus is determined by the relation Lin = /m-^d- The BLR radiation held relies on the 
radial distance of the emitting region from the BH, thns we can relate to the location of 
the emitting region via ~ tan^j. Since r^ — Sx 10“^ pc which is less than the scale of the 
BLR, the VHE photons will be absorbed by BLR photons. Thus, we must adjust the radial 
distance to avoid the absorption of VHE photons, and hnd that the radial distance 
satishes the SED modeling well. 

The modeling SEDs are shown in Fignre 3, and the parameters are listed in Table 1 
and Table 2, in which the cold proton power is estimated by assnming one cold proton per 
relativistic electron. It should be emphasized that the radiative power given in Table 2 is 
absolute, which is related to apparent one by (see Appendix in Dermer et al. 2012) 

Tabs,syn/SSC Skj Lgy^i/SSC 

Tabs,EC = 32rj^LEc/5(5D. (6) 

It is noted from Fignre 3 that the low-state SED from 2008 Angust 4 to September 12 
is well reprodnced by the SSC model, showing no hint of Compton scattering of BLR and/or 
dusty torus photons. This indicates that the location of the emitting region in the quiescent 
state is well beyond the dusty torus, at least r > Rtoms — 4.8 pc. As for both flaring states, the 
SEDs show signihcant Compton dominance, in which the parameter Aq where 

and are isotropic Compton and synchrotron Inminosity) are 29.5 and 28.5 for SED-2 and 
SED-3, respectively. While the ratios of isotropic SSC and synchrotron radiative Inminosity 
are 11.4 and 12.4 during the flaring periods. The SED is comprised by Compton scattering of 
BLR and dusty torus radiation helds, which have comparable flux levels. Table 1 and Table 
2 show that both flaring states have almost the same parameters, implying that they prodnce 
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radiation under roughly similar circumstance. 

3.3. Comoving energy density of the radiation fields 


Because the locations of the emitting regions of 4C +21.35 are within the BLR clouds 
during the flaring states, the cooling of relativistic electrons becomes more complicated. Here, 
we explore which cooling process dominates particle cooling in the emitting region. 

The energy-loss rate of electrons in Thomson limit depends on the energy density in the 
comoving frame: 


IL 


ext 


= j) dQ'J de'{ —)\(e*,r2*). 


(7) 


Here, the photon energy e* and angle 0* = cos“^p* in the stationary frame of the accretion 
disk are related to the comoving photon energy e' and direction cosine p' through the relations 
e* = rje'(l + + /^r)/(l + and 0* = f. Therefore, the comoving energy 

density of accretion disk is 

3 Re,Ld /’^max dfi' ip{r) 


Wdiscl’^) = 


16 ncgTf + 


( 8 ) 


where cfir = c/y — 1 is the velocity of the moving plasma. 

In contrast to the emission from accretion disk seen in the comoving frame, the emission 
from BLR will be Doppler-boosted, the total comoving energy density for emission lines is given 
by 


■^^BLR.lineV J ~ 


27rV'^ ™ 


N, 


e*,i 




d6^, sin 0* 


X Jiine(r,6>*)(l -/3rcos6>*)^ (9) 

where m is the number of the emission lines. The total comoving energy density of the contin¬ 


uum IS 


'^^BLR.contvJ ” 




dR sin d* 


X /cont(r,6'*)(l -/3rcos6»*)^ (10) 

At last, the total comoving energy density of the IR dusty torus can be easily written as 
- . X /iR-hd r 1 
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IR 


r = 


247rRiVrn(i + ^^^/ 


(l+/3rf*U)=J' 

The energy densities of the external radiation helds in comoving frame of the relativistic 
jet as a function of radial distance r from the BH are shown in Figure 4. Close to the BH, at 
r ^ 0.05 pc, the radiation field from accretion disk is the main target source for ERC process. 
As the distance increases, the energy density of the accretion disk rapidly decreases. At larger 













distance the BLR starts to dominate the cooling of electrons, and up to outer radius of the 
BLR, i.e., r = rsLR.out — 0.26 pc. At r > rBLR,out the radiation held from dusty torus dominates 
the electron cooling. Thus the dominant radiation held responsible for the electron cooling 
depends on the distance of the emitting region from the BH. For 4C +21.35 the emitting 
regions are within the BLR clouds, the BLR photons will play a dominant role than the dusty 
torus radiation. 

In addition to external radiation helds, the synchrotron radiation also cools elec¬ 
trons. The comoving energy density of the synchrotron radiation is related to the appar¬ 
ent isotropic synchrotron luminosity via ~ ~ 5.33erg cm“^, larger than 

Ublr — 2.63erg cm“^, which is the energy density of the BLR at r ~ 0.25 pc where the haring 
states take place. This indicates that the synchrotron radiation dominates electron cooling even 
though the BLR radiation is Doppler-boosted in the comoving frame. 

4. Discussion 

The discovery of the VHE emission between ~ 70 and 400 GeV with hux-doubling time 
about 10 minutes from 4C +21.35 is very important to limit the location of the emitting region 
and the structure of the relativistic jet. The 10 minutes variability constrains the size of the 
emitting region to be ~ 5 x 10“^[(5d/30] pc. If the size is related to the cross section of 
relativistic jet with a half-opening angle ~ 0.1, the location of the emitting region well within 
the BLR cavity, implying that the jet of the 4C +21.35 is stratihed, in which the compact 
emitting region is only a portion of the spine. 

Observational investigations based on a substantial sample of blazars using VLBI imag¬ 
ing suggest that the bulk Lorentz factor and the half-opening angle of the relativistic jet satisfy 
~ 0.1 —0.2 (Pushkarev et ah 2009; Clausen-Brown et al. 2013). Theoretically, numerical 
simulations of acceleration and collimation of external-pressure-supported relativistic jets also 
hnd that rj6'j ^ 1 (Komissarov et ah 2009; Tchekhovskoy et al. 2010). Therefore 6*j ~ 0.2/rj is 
adopted in our SED modeling, while Tj ~20 from apparent superluminal measurement (Jorstad 
et al. 2001), this results in 9-^ ~ 0.6. Our modeling hnally determines the location of the emit¬ 
ting region at rt = 8i7[,/tan6'j. Actually, we can obtain rt = 4i?(,/tan6*j through the relation 
6 *j ~ O.l/Tj. This hne-turning will increase the value of the Doppler factor and slightly decrease 
the normalization factor of the electron distribution, but has no signihcant influence on other 
parameters. Thus we can use the commonly relation 9-^ ~ (0.1 — 0.2)/rj to constrain the half¬ 
opening angle of the jet of 4C +21.35, which is within the range 0.29° ^ 9j ^ 0.6°, implying 
that the emitting regions locate around the outer radius of the BLR. 

In the optical/ultraviolet band, the clear hard spectrum and the absence of the long¬ 
term flux variations favor the dominance of a thermal component originating from the accretion 
disk. We thus model these data points with a multi-colored Shakura-Sunyaev disk, assuming 
A^bh = 6 X 10® Mq and Edisk.max = 10 eV. The modeling SED to optical/ultraviolet data is shown 
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in Figure 3 with luminosity = 4.6 x 10^® erg s“^. We note that F'disk.max = 10 eV only gives 
a good ht for SED-3, we must adjust Edisk.max to ht SED-1 and SED-2. Figure 5 shows an 
exclusive modeling to optical/ultraviolet data by varying the peak energy of the disk from 
10 eV to 16 eV. From the statistical and systematic uncertainties, the spectrum of disk could 
be described by peak energy Edisk.max = 14 eV. 

Generally, the maximum temperature of the accretion disk, corresponding to the maxi¬ 
mum photon energy, is assumed to be at the innermost radius. We can calculate the innermost 
radius with the peak energy during SED modeling as 


pmin 

-^disk 


Rp 


-13.3 


^d,46 


-, 1/3 


-M| (Edisk,max/10eV)4 r/d/(l/12) ’ 

where is gravitational radius and r/d = 1/12 is the efhciency of the gravitational potential 
energy of accreting matter transforming to the radiation. Using our modeling results that 
Ld = 4.6 X 10"^® erg s“^ and Edisk.max = 14 eV, the innermost radius is 0.22i?g, less than the 
normal innermost stable radius 6i?g of the 4C -1-21.35, this results apparently is unrealistic. 
Even if we take Edisk max = 12 eV, which corresponds to 1.86 i?g, this results remain problematic! 
If Edisk max = 10 eV, as used in our all-multiwavelength SED modeling, the innermost radius is 
~ 6.7Rg, which agrees well with the innermost stable radius of a standard Shakura-Sunyaev 
disk. This argument provides a further support for the BH mass of 4C -1-21.35 with Mbh = 
6 X 10® Mq. On the other hand, this result also indicates that the ultraviolet data points shown 
in SED-1 and SED-2 are problematic due to some certain reasons, the use of these data points 
will lead to unrealistic innermost radius of the accretion disk. 


5. Conclusions 


4C -1-21.35 was identihed as a VHE y-ray emitter by MAGIC on 2010 June 17, where 
the source was in a period of high y-ray activity detected by Fermi-LAT. The de-absorbed VHE 
spectrum smoothly connects to MeV-GeV spectrum from Fermi-LAT at ~ 3 GeV. In the paper, 
we study the location of the emitting regions using three quasi-simultaneous SEDs, including 
a quiescent state of 2008 September 12 and two flaring states of 2010 April 29 and 2010 June 
17. We preform the SED modeling based on homogeneous one-zone leptonic model, in which a 
flat BLR and a spherically symmetric dusty torus contribute external photon helds. We obtain 
the following results: 

(i) The relativistic jet of the 4C -1-21.35 is conical with a half-opening angle within the 
range of 6*j — 0.29° — 0.6°. 

(ii) The emitting regions locate within the BLR clouds and around the outer radius 
of the BLR during flaring states, while the emitting region locates beyond the dusty torus, 
r > Rtorus — 4.8 pc, for quiescent state. 

(iii) The quiescent high-energy emission is dominated by SSC process. The high-energy 
emission during flaring periods is dominated by Compton scattering of the BLR and dusty 
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torus radiation fields, which have comparable flux levels, the values of the Compton-dominance 
parameter are 29.5 and 28.5 for SED-2 and SED-3, respectively. 

(iv) The SED modeling to optical/ultraviolet data supports that the BH mass of 4C 
+21.35 is Mbh = 6 X 10® M©. 
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Fig. 1. Pair-production opacity r versus the aperture angle a of the BLR for different 7 -ray energies from 
Ery = 50 GeV upward to = 400 GeV. The location of the 7 -ray emitting region is fixed at rb = 5 x 10“^ pc 
from BH, which corresponds roughly to 10 minutes VHE 7 -ray variability. The accretion disk luminosity 
is Ld = 5 X 10^® erg s“^. The horizonal dashed line corresponds to r = 1. 



a [deg] 


Fig. 2. The critical position Rc of the 7 -ray emitting region versus the aperture angle a of the BLR for 
different 7 -ray energies from E^ = 50 GeV upward to Ej = 400 GeV. The shaded region represents the one 
surrounded by inner radius rBLR.in and outer radius r-BLR.out of the BLR. The accretion disk luminosity 
is Ld = 5 X 10^® erg s“^. The horizonal thick line represents the scale of the dusty torus. 
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Fig. 3. Spectral energy distributions of 4C +21.35 in three epochs: 2008 August 4 to 2009 September 
12 (dark-gray diamonds), 2010 April 29 (blue squares) and 2010 June 17 (red circles). Dashed lines 
indicate synchrotron, the dusty torus, the accretion disk and SSC emission components, as well as the 
Compton-scattered the dusty torus and the BLR emissions, respectively. The thicker solid line is the sum 
of all the components. The MAGIC data have been corrected for the EBL absorption using the model of 
Finke et al. (2010). 


Table 2. The comoving energy density of the electrons (f/g) and magnetic field (f/g); The jet power in the form of bulk 
motion of electron (Pe), proton (Pp) and Poynting flux (Pb); Absolute luminosity of synchrotron (Lgyn), SSC (Tssc) as 
well as Compton-scattered dusty torus (P^) and BLR radiation fields; The total absolute radiation luminosity (Pr) 

and the total jet power (P*°*). Energy density is in units of erg cm“®. Power and luminosity are in units of erg s“^. In 
calculation to the power related to protons, assuming one cold proton per relativistic electron. 
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Fig. 4. Comoving energy density of accretion disk, flat BLR and IR dusty torus radiation fields 
versus the radial distance r from the BH. The shaded region is surrounded by inner radius reLR.in and 
outer radius rBLR.out of the BLR. The vertical dashed line within the shaded region represents the location 
of the emission regions during both flaring states. 
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Fig. 5. Similar to Figure 3, but zoomed in optical and ultraviolet portion of three SEDs, which in our 
model originates mainly from a Shakura-Sunyaev disk emission. Model disk emission for several peak 
energies Edisk.max between 10 eV and 16 eV (with increment of unity) from left to right are shown. 
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